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Abstract. We revisit the published analyses of ROSAT X-ray observations of the star forming regions NGC 2264 
and Chamaeleon I (~ 3 and 5 Myr old respectively) in the light of newly published optical data. At odds 
with previous results on Chamaeleon I members, we find that low mass stars in both regions have near-saturated 
emission levels. Similarly to what previously found in the Orion Nebula Cluster, Weak Line T-Tauri Stars in 
NGC 2264 and in the Chamaeleon I cloud have higher X-ray activity levels respect to Classical T Tauri Stars, 
arguing in favor of a role of the disk and/or accretion in determining X-ray emission. 
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1. Introduction 

The influence of accretion disks surrounding young PMS 
stars on their observed X-ray activity levels is presently 
debated. The topic has been investigated many times 
through imaging X-ray observations of star forming re- 
gion but contradictory results are reported. Mentioning 



just a few recent examples, Feigelson et al. (2002) analyze 



Chandra ACIS-I data finding no indication that the pres- 
ence of an accretion disk modifies activity levels of Orion 
Nebula Cluster (ONC) stars. The same negative result, 
although so mewhat controversial, i s repo r ted fo r IC 348 
members by Preibisch fc Zinnecker ( 2001 , 2002 ), al so us- 



ing Chandra AClS-l observations; by Lawson et al. ( 1996 ) 
for the Chamaeleon I c loud using ROSAT PSPC data; by 
Flaccomio et~aLl (|2000| ) for N GC 2264 using the ROSAT 
HRI; by Grosso et al.| (200C) for p Ophiuchi again with 
the ROSAT HRI; by [Getman et al.| ( |2002| ) for NGC 1333 
(ACIS-I). 

On the other hand, Classical T-Tauri Stars (CTTS) 
belonging to the Taurus- Aurigae association are found 
to be sub-luminous in the X-ray band respect to Weak 
Lined T-T auri Stars (WTTS) b y both [Neuhauser et al 
(| 1995|) and [Stelzer fc Neuhausci| ( |2001| ). [Flaccomio et aT 
(2002b), using Chandra HRC-I data, report a similar re- 
sult, with high stati stical confidence and at odds with 
Feigelson et al. ( 2002 ), for the rich ONC population. Other 
indications of a difference between CTTS and WTTS have 
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been found in X- ray band var i ability characteristics and 
spectra. Namely, Stelzer et al. (2000) in Taurus- Aurigae, 
Flaccomio et al.| (^000|) in NGC 2264 and Flaccomio et 
al. (in preparation) in the ONC, all find that CTTS are 
more variable than WTTS. Some studies have also indi- 
cated that CTTS may have d ifferent X-ray sp e ctral char- 
acteristics respect to WTTS: Tsujimoto et al. ( ^002 ) find 
that the mean kT for CTTS is about 3 keV, compared to 
1.2 for WTTS. Such a large kT difference may in part 
be due to a selection effect: in the X-ray selected sample 
of Tsujimoto et al. (2002) class II sources (CTTS) are sig- 
nificantly more absorbed respect to class III-MS sources 
(WTTS) and it is therefore possible that only the hard- 
est CTTS have been observed. Other contrasting indica- 
tions have been also presented: Kastner et al, ( ^002 ), using 
high resolution X-ray spectra of the lOMyr old CTTS TW 
Hydrae, derive a differential emission measure distribu- 
tion peaking at ~ 0.3 keV and propose that the emission 
mechanism is related to matter accretion. No systematic 
di fference in kT between CTTS and WTTS is observed 
by Preibisch & Zinnecker ( 2002| ) in IC 348 members. 



Are these contradictory results due to real differences 
between different star forming regions or to the differ- 
ent approaches used in analyzing and interpreting data? 
We will touch upon four important points that can affect 
the result: 1) accounting for the mass/Lboi dependence 
of PMS activity; 2) choosing a relevant accretion/disk in- 
dicator; 3) avoiding selection effects in the reference stel- 
lar sample; 4) converting observed X-ray photon detection 
rates to X-ray luminosities. 
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Mass/Lbol dependence. It is possible that the failure to 
detect a difference in activity levels between stars with 
different circumstellar/accretion properties is due to the 
fact that the activity levels are also influenced by other 
stellar characteristics, and the various contributions have 
not been disentangled. In particular a dependence of mean 
Lx on stellar mass (or bolometric luminosity) has been 
widely found for PMS stellar groups. Most of the studies 
mentioned above compare the X-ray Luminosity Function 
(XLF) of CTTS with that of WTTS, both XLFs being 
computed from stellar samples comprising a wide range 
of masses. Such a procedure tends to hide possible differ- 
ences because: 1) the presence and/or magnitude of the 
effect sought might depend on stellar mass; 2) if a different 
mass-Lx relation holds for these two classes, stars having 
the same Lxj but different mass, will contribute to both 
XLFs. A more sensible approach, in order to eliminate 
this source of confusion (see e.g. Flaccomio et al. 2002b| ), 
is to compare XLFs of stars in restricted mass ranges; this 
however requires, for meaningful statistical comparisons, 
sufficiently large samples of well characterized members. 
To the same effect, considering that the ratio Lx/^boi is, 
for low mass PMS stars (< 3Mq), on average much less de- 
pendent from mass than Lx (e.g. Flaccomio et al. 2002b| ), 
it is also sensible to compare the distributions of this lat- 
ter parameter for the two classes. As an added advantage, 
Lx/Lhoi is arguably less sensitive to interstellar extinction 
corrections (although the newly introduced variable, Lbob 
also carries non- negligible uncertainties) . Both of these ap- 
proaches were successfully followed by p^"laccomio et al 
(2002b) to establish the difference in activity levels be- 
tween accreting and non-accreting ONC members. 



Disk/accretion indicator. There is no widespread consen- 
sus on which indicator of presence of disk or of accretion 
is to be used to search for effects on activity levels. Some 
studies have correlated X-ray data with accretion indi- 
cators, such as the Ha or Ca II line emission. Others 
have instead employed circumstellar disk indicators such 
as near IR excesses (in K or L band). The relation between 
presence of disks and matter accretion phenomena is not 
yet fully understood. A statistical correlation between ac- 
cretion and disk indicators is generally observed but it 
seems clear that not all IR detected disks are associated 
with accretion and it is also possible that the presence 
of accretion is not always related to a near-IR detectable 
disk (e.g. because of a large inner disk hole that suppress 
the K band excess or because of the disk view-angle). It 
is presently unclear whether X-ray emission levels are in- 
fluenced by one (or both) of the two phenomena. It seems 
therefore reasonable to explicitly distinguish between the 
two. As a clarifying example, Preibisch & Zinnecker (2002) 
find statistically significant evidence that accreting stars 
in IC 348 have lower Lx respect to non accreting ones, but 
no evidence of a difference between stars with and with- 
out K — L color excess (a disk indicator). They consider 
these two results contradictory, maintain that the infrared 



excess gives a more realistic picture of the circumstellar 
properties of the T Tauri stars than the Ha emissior^ 
and attribute the detected difference in X-ray luminosity 
functions to selection biases. An alternative point of view 
would be that accretion and not the presence of an IR- 
detectable disk influences X-ray activity and the two are 
not simply related. We note that in other cases in which 
a difference in activity levels is reported the distinction 
was performed on the basis of accretion indicators: Ha in 
Taurus (e.g. Stelzer & Neuhauser 2001) and Ca II in the 



ONC (Flaccomio et al., 2002b) 



Parent sample selection. Ideally, a complete and not- 
contaminated sample of members of a given SFR should 
be used to investigate the matter. However, this is in 
practice hardly possible. Understanding selection biases 
is therefore crucial: in t he cas e of IC 348 discussed above, 
Preibisch fc Zinncckci ( 2002| ) suspect that strong, easily 
detectable Ha emission may favor the inclusion in the ref- 
erence stellar sample of optically (and X-ray) faint accret- 
ing stars and therefore artificially depress the mean Lx of 
CTTS. This, i.e. selecting members on the basis of their 
circumstellar/accretion properties, or in any other way 
that favors the selection of faint CTTS over that of WTTS, 
is indeed the main risk to be avoided, or accounted for, in 
such a study. Taking the approach of dividing the whole 
sample in narrow mass ranges (see above), this problem, 
usually worse for low mass, low Lx stars, is reduced for the 
higher mass ranges. Other member selection methods, are 
not likely to result in spurious results: selection through 
sensitivity limited X-ray observations, for example, will 
sample to the same minimum Lx both CTTS and WTTS 
(assuming similar X-ray spectra and absorptions), so that 
if the two underlying (i.e. complete) XLFs do not differ, 
the detection fraction of both classes will be the same and 
the two distributions of observed Lx will not differ ei- 
ther. Inclusion in the reference sample of contaminating 
non-members, usually low Lx stars, will depress the mean 
activity levels inferred for non accreting/disk-surrounded 
stars, therefore going in the direction of producing the 
opposite result respect to the observed one. 

X-ray count-rate to Lx conversion. X-ray telescopes detect 
photons in wide energy ranges. The conversion between 
detected count-rate and X-ray luminosity depends there- 
fore on the incoming spectra. Given the low statistic of 
most X-ray sources with present day instruments and/or 



The reason Preibisch & Zinnecker (2002) do not trust H, 



as an indicator are its time variability and the fact that part 
of the Ha flux can be of chromospheric origin. However vari- 
ability should if anything tend to lower the significance of the 
difference between CTTS and WTTS and the chromospheric 
origin of part of the Ha flux should at most produce an ef- 
fect opposite to that observed: the sample of strong Ha stars, 
interpreted as accreting stars, will be indeed contaminated by 
chromospherically active stars, which are also likely to have 
high coronal activity and thus Lx- 
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the lack of spectral resolution of some X-ray detectors, 
it is necessary to assume a source spectrum.^ For coro- 
nal sources this usually reduces to assuming the kT of 
a thermal emission spectrum and the hydrogen column 
density, iVjj, of the absorbing material between the source 



and the observer. Differences in the way these two pa- 



rameters are estimated can lead to significantly different 
conversion factors and therefore affect the result of our 
search for a difference between CTTS and WTTS. On 
one hand, systematic differences in the X-ray spectra {kT 
and iVn) of CTTS and WTTS might result, if not prop- 
erly accounted for, in spurious results regarding the dif- 
ferent X-ray luminosity of the two classes.^ On the other 
hand, significant random errors in the conversion could 
easily wash out an existing correlation between Lx and 



accret ion/ disk indicators . In th e works of Feigelson et al 
( 2002 ) and Gctman et al. ( ^002 ), for example, the evidence 



presented against a difference of WTTS and CTTS is ob- 
tained using detect or plane (non a bsorption- corrected) X- 
ra y fluxe s. Likewise Lawson et al. ( 1996| ) and Flaccomio et 
al.^20(i^), also obtaining a negative result, assume a single 
count-rate to flux conversion factor and therefore neglect 



any difference in absorption between sources. Flaccomio 
et al. ( iGG2a ^) on the other hand, although assume a sin- 



gle kT for all sources, correct for individual absorption 
values (iVn oc A y), finding a positive result, ^telzer fc 
Neuha^se^ ( ^001| ) also correct for absorption (through a 
low energy hardness ratio) and find a positive result. Is 
it possible that a peculiarity in the A^h vs. Ay relation 
or in the intrinsic spectra of CTTS and WTTS results 
in artificially lowering the luminosities derived for CTTS 
respect to those derived for WTTS? Regarding the rela- 
tion between X-ray and optical extinction, recent stud- 
ies correlating A'h, derived from X-ray medium resolution 
spectra, and optically derived Ay confirm the relation be- 



Imanishi et al.. 


2001; 


Flaccomio et al, , 


2002a; 


Feigelson et 


al., 


20( 


12; 


Kohno et al. 


2002). It is however possible, al- 



though presently not still fully established, that CTTS 



have harder X-ray spectra respect to WTTS ( Tsujimoto 
et al 



2002) 



We may wonder how a different kT would 
affect the luminosities we derive. Figure^ shows the value 
of the conversion factors for the ROSAT HRI, as a func- 
tion of kT and A'h- Qualitatively similar plots are ob- 
tained for the ROSAT'PS^G and for the Chandra HRC-I, 
the two other instruments used for the observations dis- 
cussed later in this paper. We observe that, for a given 
source A^h, the difference in kT, if eventually confirmed, 
will indeed go in the direction of decreasing the inferred 
Lx of lower kT sources respect to high kT ones, thus po- 
tentially accounting for part of the observed differences 



^ Even in the case of detectors with intrinsic spectral resolu- 
tion, the low source counts practically prevent in many cases 
the determination from the data of a reliable spectral model. 

^ Differences in the intrinsic spectra would however be inter- 
esting "per se" for the understanding of the physical mecha- 
nism that determines X-ray emission. 



between CTTS and WTTS. However it is also clear that 
for typical Ay ~ 0.5 — 4.0 the mistake committed in not 
accounting for individual source temperatures could be at 
most of the order of ~ 0.1 — 0.2 dex, smaller than the 



difference bet ween CTTS and WTTS found by [Flaccomic 
et al.| ( |2002b| ) in the ONC. We stress however that even 



if this effect were to be confirmed, therefore reducing the 
actual difference in luminosities respect to that inferred 
assuming a single fcT, the difference between the X-ray 
emission of CTTS and WTTS would be confirmed, and 
the spectral differences would provide additional clues for 
the understanding of its physical origin. 

In this paper, keeping the above four points in mind, 
we further discuss and extend the evidence for a role of 
accretion and/or disk in determining the observed X-ray 
activity level of ONC members, as already reported by 
Flaccomio et al. (2002b). In the light of newly available 
optical/IR data wc then critically reanalyze the results 

( ^000 ) and iLawson et al 



obtained by Flaccomio et al 



( 1996D , both of which concluded that stars surrounded by 



disks, in NGC 2264 and Cha I respectively, have the same 
activity levels as those that do not have a disk. Here we 
derive the opposite result. 

The structure of this paper is as follows: in Sect. || we 
discuss the new observational evidence for a difference in 
activity levels between CTTS and WTTS belonging to the 
ONC. In Sect. | and Sect. | we then discuss the cases of 
NGC 2264 and the Chamaeleon I cloud. Finally in Sect. 
H we briefly summarize our results. 

2. The ONC: an update 

Our analysi s method throughout th is paper is based on 
the work of [Flaccomio et al.| (|2002a| ,|b|) on the ONC. This 
cluster is arguably the best target available for our study 
because we have access to a rich and well characterized 
sample of member s spanning a wide range of masses. We 
refer the reader to Flaccomio et al. ( 2002a ^ for a full de- 
scription of the X-ray and optical data used here. The ex- 
tinction limited sample (Ay < 3.0) discussed in this latter 
work has little field contamination and is complete almost 
down to the lowest stellar masses. Using the Ca II line 
(A — 8542A) as an ind icator for circumstellar accretion, 
Flaccomio et al. ( 2002b ) obtained with high statistical sig- 
nificance the result that low mass stars (M < 3Mq) with 
this line in strong emission {EW < —1) have systemati- 
cally lower Lx and Lx/iboi values respect to stars with 
the line in absorption {EW > 1). 

Here we state that an analogous result is obtained, 
albeit with smaller significance (2.5 — Su), comparing the 
Lx and Lx/L]^oi distributions of stars with large and small 
near IR excess (A(/ - K) > 0.8 and A(/ - K) < 0.4 re- 
spectively). The X-ray and optical/IR data are presented 
in Flaccomio et al. (p002a ) and the A(/ — K) values are 
taken from Hillenbrand et al. (1998). Figure]^ and|3|shows 
the maximum likelihood Lx and Lx/iboi distributions for 
these two classes of stars in 6 different mass bins. The 
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surements, extinction toward sources had to be assumed 
uniform and placement of counterparts in the HR diagram 
was performed solely from optical photometric data, an 
error-prone procedure for PMS stars. Moreover the dis- 
tinction between CTTS and WTTS was not well estab- 
lished, as indications on the NIR excesses were not avail- 
able and Ha measurements were in most cases qualitative 
and non-uniform. Since that work new improved optical 



data have been recently published by RebuU et al. (2002) 



Fig. 1. Contour plot of the logarithm of the conversion 
factor {CF) from ROSATHYd count-rates to unabsorbed 
flux, in the 0.1-4.0keV spectral band, for coronal thermal 
sources (Raymond-Smith emission model) as a function of 
source temperature, fcT, and absorption, A^h- Lines refer 
to loci in which Log{CF) is constant and equal to the 
value reported in the corresponding label plus -10.07, i.e. 
the Log{CF) corresponding to typical parameter values: 
kT = 2.16, iVH = 2 • 10^1 (marked in the Figure by a tri- 
angle). Line spacing is 0.1 dcx in Log{CF). The ordinate 
axis on the right gives the standard conversion between 
A'h and optical absorption. Ay. iVn = 2 • lO^^^v- 

range of confidence with which we can exclude that the 
two distributions are randomly extracted from the same 
parent population, according to the tests in the ASURV 
package (Feigelson & Nelson, 1985| ), is given inside each 
panel. 

We also note that very similar results are obtained, 
both with the Ca II line and with A(/ — K) as a discrim- 
inant, if only X-ray detected members are considered in 
the distribution functions. As reminded in the introduc- 
tion, these latter results exclude that the difference in the 
distributions is due to the preferential selection of faint 
CTTS (which is anyway not expected given that the sam- 
ple is noi selected from either accretion or disk indicators). 

3. NGC 2264 

Stellar activity in NGC 2264 has been studied most re- 
cently by Flaccomio et al. ( 2000D through 6 different 
R0SATY{J{1 observations covering, in two different point- 
ings, a large fraction of the star forming cloud. One hun- 
dred sixty nine distinct sources were detected, ~ 95% of 
which are estimated to be associated with members of the 
association. One of the main problem at the time of this 
work was the lack of good optical characterization of the 
members, so that, for example, lacking individual mea- 



We therefore updated the previous analysis according to 
the general principles stated in the introduction. 

3.1. The reference sample 



We a d opted optical data from tables 1 and 3 of [Rebul] 
let al.| (|2002| ). Out of the full hst of 687 photometrically 
selected candidate members (i.e. the in cloud, in locus 



population defined by RcbuU et al. 2002) we selected the 
202 stars for which reliable spectral types and extinction 
{Ay) estimates (through R and / photometry -I- spectral 
types) were available. This latt er is a subset of th e full 
spectroscopic sample studied by RebuU et al. (2002), ini- 



tially selected primarily from a list of I-band variable stars, 
with the addition of previously known candidate members 
based on their X-ray or Ha emission or on their proper 
motion. Our reference sample is therefore the intersection 
(logical and) of the photometrically selected member sam- 
ple and of the spectroscopic sample. While the former is 
arguably free from selection biases in favor of faint CTTS, 
the degree of representativeness of the latter in this re- 
spect is less clear: Ha is however only a secondary selec- 
tion criterion and I-band variability (periodic in > 50% of 
the stars), although maybe more frequent in CTTS, does 
not obviously favor the inclusion of optically fain t stars . 
Moreover the disk (CTTS) fractions [RcbuU ct all ( |2002|) 
derive for the spectroscopic and the photometric sample 
(cf. their Table 6) are remarkably similar, suggesting that 
the former is not strongly biased toward CTTS. 

Contamination of our reference sample from field stars 
may on the other hand be non negligible: according to pre- 
liminary proper motion data RebuU et al. ( ^002 ) report 
that ^ 50% of their photometric candidates are actually 
non-members. The spectroscopic sample, selected on the 
basis of PMS stellar characteristics, is expected to be less 
contaminated, although an estimate based on proper mo- 
tion data is not provided. We recall (cf. Sect. |l|) that field 
star contamination is expected to artificially lower the ac- 
tivity levels of WTTS. 

We place stars in our reference sample in the HR 
diagram. Effective temperatures and bolometric correc- 
tions are estimated from spectral types and Kenyon & 
Hartmann ( 1995 ) conversions]^ We then evaluate bolo- 



metric luminosities from I band magnitudes. Out of the 



^ The effective temperature for spectral type A6 was modi- 
fied from 8350 K to 8050. The former value seems to be at odds 
with the spectral type-Teff relation; we suspect a typographical 
mistake. 




Fig. 2. X-ray luminosity functions for stars with high- and fow-NIR excess (soHd and dashed hnes, respectively) in the 
Orion Nebula Cluster. Panels refer to different mass ranges, as indicated by legends. Also reported are the numbers 
of detected (d) and undetected stars (u) used for XLFs of high- and low-accretion subsamples and the cr-equivalent 
significance range for the difference between the two distributions (see text). 



202 spectr ally characterized candidate members, 193 fall 
within the 3iess et al. (200C) evolutionary model grid and 
have therefore been assigned a mass and an age. 

The reference sample used for our following analysis 
comprises the 178 stars, out of these 193 candidate mem- 
bers characterized in terms of mass and age, that fall in 
the field of view of th e X-ray observations described by 
(|2000|). 
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3.2. X-ray data 



We matched the photometric catalog of RebuU et al 



(2002), out of which our reference member list is drawn 
with the list of 169 X-ray sources published by Flaccomio 
( 20001) . The identifications were carried out as de 



et al 



scribed in Flaccomio et al. ( 200C| ), i.e. assuming as identi- 



fication radii the the off-axis dependent X-ray source posi- 
tion error summed in quadrature to 1", i.e. a conservative 
estimate of the optical position error. Before performing 



the final identification we first registered the coordinate 
systems of the optical and X-ray lists by comparing the 
positions of 125 uniquely identified pairs (i?Ax — -R^opt = 
— 1.7", Decx — Decapt = 1.9"). Sixty seven stars in our 
reference sample were identified with an X-ray source, 56 
of which uniquely, while the the remaining 11 fell in the 
identification circle of an X-ray source along with other 
objects in the photometric catalog. To each of the 67 can- 
didate members with X-ray counterparts we then assigned 
a Maximum Likelihood (ML) X-ray count rate: these are 
values computed by Flaccomio et al. (2000) in order to 
define a mean source brightness among 6 different ob- 
servations. Due to source variability, five of these mean 
values are actually upper limits. The count rates of the 
11 ambiguous identifications were also treated as upper 



limits and upper limits, computed as in Flaccomio et al 
( pooo| ), were also assigned to 111 X-ray undetected candi- 



date members lying within the FOV of the HRI observa- 
tions. 
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Fig. 3. Same as Figure for Ly^/ L 



bol- 



Finally we converted count-rates, measured and up- 
per limits, to X-ray luminosities in the 0.1-4.0 keV band|. 
We assumed a thermal emission spectrum with kT = 
2.16 keV, close t o recent estimates for PMS stars (e.g . 



Flaccomio et al. , 2002b ; [manishi et al. , ^001 ; Getman et 



al. , 20C 2 ) . The hydrogen column density was assumed pro- 



portional to the optical extinction measured individually 
for each star: A'h = 2 • IQ^^Ay- The distance was assumed 
to be 760pc like in Flaccomio ct al. (2000). Figure || com- 
pares the Lx derived by Flaccomio et al. ( 2000| ) to those 
derived here from the same count-rates. Our new estimates 
are ~ 0.3 dex higher respect to the old ones, with the main 
differences due to the assumed kT (2.16 vs. 0.75 keV) and 
the increased absorption (mean Ay ~ 0.45 vs. a constant 
Ay — 0.19) and a small difference, < 0.1 dex, due to the 
different spectral band in which Lx is computed. 



Conversion factors to un-absorbed flux are com- 
puted using PIMMS (Portable, Interactive, Multi- 
Mission Simulator) , version 3.0, available on-line at: 



3.3. Activity vs. circumstellar environment 

We will now try to establish whether the activity levels of 
our sample of NGC 2264 candidate members depend on 
circumstellar properties, evidenced either by accretion or 
by IR disk indicators. We will use two indicators provided 
by RebuU et al. ( 2002| ): Ha equivalent width, EWna, and 
excess in the H — K near IR color, A(i? — K). In the 
former case we will consider as accreting those stars with 
EWyio > 5 and non- accreting those with EWna < 5^. In 
the latter case we take as threshold A{H — K) = 0.15 with 
stars showing larger excesses considered as surrounded by 
circumstellar disks. 

Figure 1^ shows the scatter plots between Lx and stellar 
mass and ix/iboi and mass, for the accreting and non- 
accreting stars. Similarly to what seen in other SFRs, Lx 



http:/ /asc. harvard.edu/tooIkit/pimms.jsp 



^ We choose 5 instead of the more customary 10 as EW 
threshold because we so obtain a better distinction between 
activity indicators of the two classes. As noted in the introduc- 
tion lowering the threshold could, if anything, result in a reduc- 
tion of the inferred difference in activity levels; the most active 
WTTS might indeed have strong chromospheric Ha emission 
and be mistaken for CTTS. 
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Fig. 4. Comparis o n of X-ray luminosities computed by 



Flaccomio et al 



(|200q ) for NGC 2264 stars (see text) 
and those recomputed from the same data in this work. 
No distinction is made here between detections and up- 
per limits. The solid line indicates the locus of equal val- 
ues; the dotted lines indicate the relation ix(This work) 
= 2 • Lx( Flaccomio et al., 2000). The scatter of points 



about the mean relation is due to the adoption of individ- 
ual extinction corrections. 



is seen to correlate with stellar mass, although the relation 
appears somewhat fuzzier in this case respect to, e.g., the 
ONC or the Chamaeleon I cases (see below). A systematic 
difference between the two classes is not readily apparent. 
However, given the large number of upper limits a more 
quantitative analysis is needed. Figure ^ shows XLFs and 
Lx/Lhoi distributions for the two Ha separated stellar 
classes, in two different mass ranges, marked in Fig. ^ by 
vertical lines, and for the whole sample. In all cases the 
distributions of accreting stars appear to lie below those 
for non-accreting ones. Statistical test (in the ASURV pack- 
age) confirm such differences with varying degree of con- 
fidence (results of the tests are given in the figure, along 
with the number of detection and upper limits that enter 
in the distributions; the same information is repeated in 
the first part of Table 0, where the results of other tests 
described below are also reported). Particularly significant 
is the difference in Lx for the whole sample (> 3.9(j). Such 
a difference might however result from a larger fraction of 
low mass, and low Lx, accreting stars (see Fig. |^). In the 
two narrower mass ranges, the result is however retrieved, 
with greater than 3cr confidence in the 0.5 — 1.0 A/q mass 
range. The difference is also observed in Lx/-^boi- Note 
that, as pointed out above, our reference sample may be 
significantly contaminated by field stars, and this would 
tend to depress the non-accreting stars distributions, thus 
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Fig. 5. Lx and Lx/Lboi vs. mass for NGC 2264 stars with 
high and low Ha equivalent width (black and and gray 
symbols, respectively). Filled circles represent detections, 
down-pointing arrows upper- limits. 

lowering the significance of the result. If we repeat the 
same analysis including only stars confirmed as members 
by their IR excess (A(/ - K) > 0.3 or A{H - K) > 0.15), 
Ha emission [EW > 5) and X-ray detection, we indeed 
find even more significant differences. Particularly inter- 
esting are the differences in Lx/Lboi, because they are less 
likely to be influenced by selection effects, a concern for 
this latter restricted sample. Table |l| reports the results of 
these tests. 

Table |l| also reports the results of the comparisons 
between the stars with and without near IR excess, for 
the same mass ranges and the two stellar samples de- 
scribed above. The same results is retrieved: stars showing 
a A{H — K) excess, indicating the presence of a disk, have 
lower activity levels respect to the complementary sample. 

4. Chamaeleon I 

Our main source of data regarding the Chamaeleon I as- 
sociation is the work by Lawson et al. ( 1996 ): we use their 
member list (117 stars in their Table Bl) and their esti- 
mates of Lboi and Tcs (available for 78 stars. Table 6). We 
adopt here a distance to the Chamaeleon I cloud of 160pc 



(Whittet et al., 1997; Wichmann et al 



than the distance assumed by Lawson et al 



1998D , 20pc larger 
(|1996D in es- 



timating bolometric luminosities. We therefore increased 
the Lboi values accordingly. Masses of 71 candidate mem- 
bers were derived from placement in the HR diagram and 
interpolation of Siess et al. ( 2000D evolutionary tracks 



The selection of candidate members in Lawson et al 
( 1996D is performed mainly on the basis of either their Ha 
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Table 1. Significance of difference between CTTS and WTTS in NGC 2264 



iVicLoS 1 i Kl Q J 


iVd(C) Nu{C) 






Q-j nn ( T ,-\r 1 


u Ly iL.yiJX. / -t^bol ) 


EW{H^) - 


Whole sample 












0.5-1.0 


5 


16 


18 


14 


3.0/3.4 


2.5/2.8 


1.0-2.0 


5 


14 


13 


10 


1.4/1.9 


1.1/2.1 


All 


18 


79 


33 


44 


>3.9 


1.5/2.1 


EW{H^) - 


Confirmed members 










0.5-1.0 


5 


16 


18 


2 


>3.9 


3.2/>3.9 


1.0-2.0 


5 


14 


13 


1 


2.1/2.7 


1.9/3.1 


All 


18 


79 


33 


7 


>3.9 


3.3/>3.9 


A{H~K) 


- Whole samph 












0.5-1.0 


1 


9 


18 


20 


1.3/1.8 


2.2/3.2 


1.0-2.0 


1 


10 


14 


11 


1.9/2.2 


2.1/2.5 


All 


4 


30 


40 


71 


2.2/2.4 


2.6/3.2 


A{H -K) 


- Confirmed members 










0.5-1.0 


1 


9 


18 


9 


1.8/2.3 


2.6/3.5 


1.0-2.0 


1 


10 


14 


2 


2.4/2.8 


2.9/3.4 


All 


4 


30 


40 


44 


2.8/3.0 


3.4/>3.9 



Note - Description of columns: (1) Mass range; (2) Number of detected CTTS; (3) Number of CTTS with upper limits; (4) 
Number of detected WTTS; (5) Number of WTTS with upper limits; (6) Range of significance (expressed in a equivalent) 
for the difference between the Lx distributions of CTTS and WTTS according to the tests in ASURV. Results greater than 
3cr are in boldface; (7) like (6), but for -Lx/iboi 



or X-ray emission. The danger of preferentially selecting 
faint s tars (both optically and in X-rays) with strong Hg 



emissic n is therefore present. However the Chamaeleon I 



association is close enough that a large fraction of interme- 
diate mass members is probably detected in the ROSAT 
PSPC X-ray observations. Also, as anticipated in the in- 
troduction, other than Lx we will also investigate the de- 
pendence of Lx/Lboi on circumstellar characteristics and, 
as a further test, we will also consider a fully X-ray se- 
lected sample. 



individual optical extinction values are taken from the fol- 
lowing works : pawson ct"al ( 1996 , Aj, Ta ble 3), Gauvin fc 
Strom] (|1992|, Ay, Table 2) , jwaltei] ( |1992| , Eb^v, Table 1) 



and jCambresy et al. ( 199§| , Ay, Table 1); whenever multi- 
ple estimates were available for a given star we choose one 
of the four values, the precedence order being the same 
as the order of citation given above. Aj and £'b-v were 
converted to Ay by multiplying by 3.55 and 3.1 respec- 
tively (Mathis, 199C). Figure ^ compares the new X-ray 



X-ray data were taken from Lawson et al. (1996): they 
quote X-ray luminosities (or upper limits) for members of 
the region, computed from ROSAT PSPC count rates in 
the 0.4-2.5keV spectral band (Feigelson et al, , 1993), using 



luminosities with those reported in Lawson et al. ( 1996[ ) 



and indicates the effects that contribute to the consider- 
able average discrepancy between the two estimates. First 
of all a difference of ~ 0.15 dex, indicated by the low- 
est diagonal thin line, is of unclear origin: we recomputed 
the conversion factor, in the 0.4-2.5keV band, assuming 
facto r: 1 PSPC count ks'^=3 x lO^^ergs s~^ . |Feigelson| kT = l .OkeV and iVn = 2.0 • lO^i, i.e. following [Feigelsor 



a constant count-rate to Lx (in the same band) c onversion 



et al. ( 1993) find that this conversion factor corresponds 
to assuming a plasma temperature kT ^ IkeV and an 
absorption by a hydrogen column, A'h, corresponding to 
Ay ~ 1. 

In order to account for differential extinction (i.e. the 
fact that star are subject to different extinctions) and to 
uniform our assumptions to the ONC and NGC 2264 stud- 
ies, we re-estimated X-ray luminosities, in our standard 
0.1-4.0 keV band. We started from PSPC count rates in 



et al. ( 1995 ), and derived a larger conversion factor, by 



the 0.4-2.5 keV band, i.e. from the Lx reported in Lawson 
et al. (|1996| ) divided by the above mentioned conversion 
factor. We then multiplied these count-rates by conver- 
sion factors between PSPC count-rates (in the 0.4-2.5 keV 
band) and luminosities (in the 0.1-4.0 keV band), com- 
puted for a kT = 2.16 keV thermal plasma emission ab- 
sorbed by an hydrogen column A^h = 2 • 10^^ • Ay and our 
assumed distance to the association (160pc). Estimates of 



~ 0.15dex, respect to the value reported by these authors. 
The other light lines show the effect of having changed 
the assumed cluster distance, the chosen spectral band, 
the plasma temperature, and the average source extinc- 
tion. The combined effects of these changes results in our 
X-ray luminosities being on average ^ 5 (0.7 dex) times 
larger than the ones formerly derived. 

4.1. Activity vs. circumstellar environment 

We adopt the distinction between CTTS and WTTS pre- 



sented by Lawson et al. (1996, Table Bl), excluding from 



our analysis 4 stars with uncertain classification, out of 
our 71 with mass estimates. The distinction is based on 
Ha emission. Our final sample comprises 28 CTTS and 39 
WTTS. 
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Fig- 7- 



Lawson ct al 



Comparison of X-ray luminosities reported by 
( 1996 ) for Chamaeleon I stars and those re- 
computed from the same data in this work. No distinction 
is made here between detections and upper limits. The 
bottom solid line indicates the locus of equal values; the 
light lines indicate the effect, on the X-ray luminosities, 
of: recomputing the conversion factor assuming kT = 1.0 
and A'h = 2 • 10^^ (see text), changing the assumed dis- 
tance, source plasma temperature, band in which Lx is 
computed and average extinction. The scatter of points 
about the highest light line is due to the adoption of indi- 
vidual extinction corrections. 

Figure || shows, with different symbols for CTTS and 
WTTS, the scatter plots of Lx and ix/-^boi with mass. 
Disregarding for the moment the difference between CTTS 
and WTTS, a trend of increasing Lx with increasing 

( |1996D and also 



mass, already noted by Lawson et al 



seen in other star forming regions, can be clearly observed. 
Lx/Lho\ seems to be close to the saturation level (10~^) 
at all masses. We note that Lawson et al. (1996), on the 
basis of their lower X-ray luminosities had excluded that 
coronal activity in Chamaeleon I members was saturated, 
contrary to what reported for other star forming regions. 
Our re-analysis of the same data shows that this result 
can be attributed in large part to the assumptions made 
in the conversion between count-rates and X-ray luminosi- 
ties and to the choice a non standard X-ray spectral band 
for the calculation of Lx- 

Figure ^ shows the Lx and ix/^boi distribution func- 
tions, separately for CTTS and WTTS, in the same two 
mass ranges investigated in NGC 2264 and for the whole 
sample. First of all we note that there is little difference 
(at the ~ Icr level) between the two XLFs referring to the 
whole population. This is indeed the same result reported 
by Lawson et al. (1996). However a look at Fig. || shows 



Fig. 8. Lx and Lx/Lboi vs. mass for CTTS and WTTS 
(black and and gray symbols, respectively) belonging to 
the Chamaeleon I region. Filled circles represent detec- 
tions, down-pointing arrows upper- limits. 

that this might be due to the inclusion of stars over an am- 
ple range of masses. If we indeed consider only stars in the 
0.5 — 1.0 A/q range CTTS appear to be underluminous re- 
spect to WTTS at the ^ Scr level, both in absolute terms 
and respect to their bolometric luminosities. ix/-^boi is 
indeed lower (at the 2 — 2.8cr level) even if we consider 
the whole sample. We obtain similar results, although of 
somewhat lesser significance, if we only consider X-ray se- 
lected stars: for example, the significance of the difference 
in the 0.5 — I.OMq range are ^ 1.5 and ~ 2. Oct for Lx and 
ix/iboi, respectively. 

As a final note we remark that less significant results 
are obtained if the same analysis is performed with the 



values of ix reported by Lawson ct al. (1996). The scatter 
of points around the mean relations observed in Fig. ^ 
as well as in the distribution functions in Fig. ^, appear 
in this case to be larger. However the difference in the 
0.5 — I.OM0 mass range remains (at the 2.2/2.8cr level). 

5. Summary and Conclusions 

We have considered three SFRs: the ONC, NGC 2264 and 
Chamaeleon I. After a critical re-analysis of the optical 
and X-ray data published in recent literature, we have 
tried to answer the question of whether stars with differ- 
ent circumstellar properties have different observed X-ray 
emission, both in absolute terms and in relation to their 
bolometric luminosities. 

In all of the analyzed cases we find that CTTS are 
underluminous respect to WTTS. This result is found in 
spite of large differences in the selection of members and 



10 



Flaccomio et al.: Circumstellar accretion and activity in the PMS 



in the optical and X-ray data used. We believe that it 
indicates a difference either in the intrinsic properties of 
X-ray emission or, alternatively, in the radiation transport 
(e.g. absorption) in the proximity of the stellar system. 
However we tend to prefer the first option: a difference 
in the relation between optical and X-ray circumstellar 
extinction, for example, might explain our result, but no 
such indication has been found to date. 

When we could investigate the matter, i.e. in the cases 
of the ONC and NGC 2264, we found that the difference 
holds both when we discriminate stars on the basis of ac- 
cretion and disk presence indicators. Therefore we are not 
able to say which of these two related aspects is most rel- 
evant for the mechanism responsible for the difference. 

Having established the reality of this effect, work re- 
mains to be done to better characterize it and to iden- 
tify its physical source. Thanks to sensitive X-ray data 
from modern X-ray space-borne telescopes and high- 
throughput optical/IR instruments this goal seems well 
within the reach of near future research. 
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Fig. 6. Distributions of Lx and ix/-^boi (left and right columns) for NGC 2264 stars with high and low Ha equivalent 
width (solid and dashed lines, respectively). Each row refers to a different mass range as indicated. Legends inside 
panels as in Fig. ||. 
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Fig. 9. Distributions of Lx and Lx/-^boi (left and right columns) for CTTS and WTTS (solid and dashed lines, 
respectively) in the Chamaeleon I region. Each panel refers to a different mass range as indicated. Legends inside 
panels as in Fig. ||. 



